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ABSTRACT 


NORDA  Code  372  is  at  present  engaged  in  research  projects  for 
extending  Airborne  Electromagnetic  (AEM)  technology  from  geophysical 
prospecting  applications  to  Navy  requirements  such  as  airborne  shal¬ 
low  water  bathymetry,  airborne  ice  thickness  determination,  and  con¬ 
ductivity  profiling.  This  report  is  part  of  our  research  program 
addressing  AEM  ice  thickness  determination  and  will  demonstrate  that 
the  Navy  requirements  for  rapid  airborne  methods  for  mapping  ice 
thickness,  in  real  time,  can  be  accomplished  with  AEM  techniques.  At 
present  there  are  not  reliable  airborne  techniques  for  measuring 
directly  sea  ice  thickness,  although  impulse  radar  has  yielded  promising 
results  by  implying  sea  ice  thickness  from  ice  roughness  measurements. 

The  preparation  of  this  technical  note  was  part  of  the  work  per¬ 
formed  for  extending  Airborne  Electromagnetic  (AEM)  technology  from 
geophysical  prospecting  to  ice  thickness  measurements  under  Program 
Element  62759N.  Sponsorship  was  provided  by  Mr.  Barry  Dillon  of 
NAVAIRSYMCOM/Code  340.  The  AEM  program  is  managed  by  Mr.  Kuno  Smits , 
NORDA  Code  372.  The  work  was  performed  during  the  period  May-August 
1983. 
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Summary 


This  repsrt  examines  the  feasibility  of  using  an  air¬ 
borne  electromagnetic  induction  apparatus  for  the  mapping  of 
sea  ice  thickness.  After  a  review  of  available  commercial 
systems  and  the  electrical  properties  of  sea  ice  it  can  be 
cc  eluded  -hat  this  novel  use  of  airborne  electromagnetics 
is  indeed  very  feasible. 

In  order  to  take  full  advantage  of  this  technique  we 
suggest  a  design  for  an  airborne  system  that  is  optimal  for 
the  designated  task  of  sea  ice  mapping.  Before  going  on  to 
its  construction  however,  we  strongly  recommend  that  the 
theoretical  ideas  presented  in  this  report  be  field  tested 
to  ensure  that  all  sources  of  error  have  been  accounted  for. 
The  required  field  test  can  be  done  with  conventional  con¬ 
tractor  apparatus. 

2 .  Introduction 

Airborne  electromagnetic  induction  systems  have  r. ow 
seen  nearly  thirty  years  of  commercial  service.  Tradition¬ 
ally,  this  type  of  geophysical  exploration  equipment  is 
optimized  for  the  detection  of  discrete,  highly  conductive 
mineral  deposits  which  occasionally  prove  to  be  of  economic 
value.  More  recently,  airborne  electromagnetic  apparatus 
has  been  used  for  geological  mapping  and  the  related 
indirect  indication  of  economic  minerals.  Because  of  the 
fair  degree  of  success  that  was  reported  in  interpreting 


1 


airborne  electromagnetic  data  in  terms  of  an  underlying  geo¬ 
logical  section  (e.g.  DeMouli.*  and  Becker  1992),  it  was  only 
natural  that  one  would  enquire  as  to  the  possible  oceano¬ 
graphic  applications  of  this  type  of  apparatus.  In  this 
context,  Morrison  and  Becker  (1982)  reviewed  the  possible 
use  of  airborne  electromagnetics  for  coastal  bathymetry. 
Similarly,  this  report  examines  the  feasibility  of  another 
possible  oceanographic  application  of  airborne  electromag¬ 
netic  namely,  the  mapping  of  sea  ice  thickness. 

3  Airborne  Electromagnetic  ( AEM)  Instrumentation 

i  •  Principles 

All  airborne  electromagnetic  (AEM)  systems  basically 
consist  of  a  primary  field  source  (transmitter)  which  pro¬ 
duces  a  time  dependent  magnetic  field  and  a  detector 
(receiver)  which  senses  the  secondary  magnetic  field  pro¬ 
duced  by  the  eddy  currents  in  the  earth.  The  first  system 
of  this  kind  was  developed  some  twenty-five  years  ago  in 
Canada  for  the  purpose  of  mineral  exploration.  Since  then  a 
large  variety  of  AEM  systems  have  been  put  into  production. 
A  majority  of  these  have  now  been  abandoned  for  one  reason 
or  another.  In  particular  at  the  time  of  the  last  review  of 
the  state  of  the  art  in  this  field  (Becker  1979),  nearly  a 
dozen  different  systems  were  available  for  survey  work.  At 
present,  this  abundance  has  basically  been  reduced  to  two 

AEM  systems,  (c.f.  Figure  1)  namely  the  towed-bird,  time- 
domain  INPUT  system  which  is  used  for  about  70%  worldwide 
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survey  work  and  the  rigid-boom  helicopter  towed  frequency 
domain  system  which  accounts  for  about  25%  of  the  work.  The 


remaining  5%  of  the  work  <s  done  with  special  purpose  and 
experimental  systems. 

The  two  available  systems  differ  widely  in  scale  but 
overlap  in  the  frequency  range  covered.  Figures  2  and  3 
graphically  illustrate  this  point.  While  the  INPUT  system 
operates  in  the  time-domain  transmitting  an  interrupted, 
high-power,  half  sine  pulse,  it  can  be  shown  by  Fourier 
analysis  that  it  covers  all  the  odd  harmonics  of  its  funda¬ 
mental  frequency  (iuu  Hz)  up  to  about  5  kHz.  The  helicopter 
system  operates  in  the  frequency  domain  and  covers  a  limited 
frequency  range.  Normally,  both  systems  detect  the  horizon¬ 
tal  component  of  the  secondary  field.  On  demand,  however, 
either  system  may  be  modified  to  measure  the  vertical  com¬ 
ponent  of  the  secondary  field  or,  if  need  be,  both  com¬ 
ponents  . 

3.2  Time  Domain  Equipment 

The  INPUT  (Induced  Pulse  Transient)  system  was  invented 
some  twenty  years  ago  by  A.  P.  Barringer  (Barringer  1962) 
and  has  been  modified  a  number  of  times  since  it  was  first 
rendered  operational.  It  uses  a  high  power  1  millisecond 
pulse  to  induce  eddy  currents  in  the  ground.  The  secondary 
magnetic  field  related  to  the  decay  of  these  eddy  currents 
is  then  detected  while  the  primary  field  is  off.  The  detec¬ 
tion  of  the  weak  secondary  fields  during  the  absence  of  any 
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FIGURE  2.  AEM  SYSTEM  GEOMETRY 


strong  primary  fields  avoids  electrical  dynamic  range  prob¬ 
lems  and  allows  the  use  of  a  towed-bird  receiver.  Simi¬ 
larly,  because  the  primary  field  is  off  while  the  secondary 
field  is  detected,  one  can  tolerate  a  fair  amount  of  motion 
of  the  detector  with  respect  to  the  aircraft  without  any 
serious  degeneration  of  the  signal  quality.  The  transient 
signal  is  sampled  at  a  number  of  times  after  the  cessation 
of  the  primary  field.  Fach  sample  is  then  fed  to  a  detector 
channel  (of  which  there  usually  are  six)  and  averaged  over 
about  150  repetitions  of  the  primary  pulse.  The  average 
value  of  each  channel  is  then  displayed  on  an  analog  record 
as  a  function  of  distance  along  the  flight  path.  A  detailed 
analysis  of  this  system  was  made  by  Becker  (1967)  while 
Palacky  (1972)  provides  a  description  of  the  most  recent 
version  of  this  apparatus. 

2*1  F r equency  Domain  Equipment 

Figure  is  a  schematic  diagram  of  a  typical  hel¬ 
icopter  frequency  domain  system.  Note  the  rigidly  mounted 
transmitter  and  receiver  which  allows  the  cancellation  of 
the  primary  field  at  the  receiver  and  thus  permits  the  pre¬ 
cise  measurement  of  the  in-phase  and  quadrature  components 
of  the  secondary  field  in  units  of  parts  per  million  (ppm) 
of  the  primary  field.  A  recent  description  of  this  type  of 
system,  mounted  on  a  helicopter-towed  9m  boom,  is  given  by 
Fraser  (1979).  Usually  as  shown  in  the  sketch,  the  hel¬ 
icopter  also  tows  a  magnetometer.  Flight  path  recovery  is 
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done  with  the  aid  of  an  on  board  film  strip  camera.  The 
helicopter  also  carries  a  radar  altimeter  so  that  its  eleva¬ 
tion  with  respect  to  the  ground  is  known  precisely.  It 
should  be  noted  however  that  the  position  of  the  EM  boom  with 
respect  to  the  helicopter  is  only  known  approximately. 

The  electrical  circuit  used  in  a  rigidly  coupled  fre¬ 
quency  electromagnetic  system  is  shown  in  Figure  5.  An 
electrical  link  exists  between  the  transmitting  and  receiv¬ 
ing  elements.  It  is  used  to  cancel  the  voltage  induced  in 
the  receiver  by  the  primary  field  and  to  provide  a  phase 
reference  for  the  detection  and  measurement  of  the  secondary 
field  . 


The  technical  specifications  for  rigidly  coupled  fre¬ 
quency  domain  airborne  systems  are  given  in  Table  1 .  The 
helicopter-towed  apparatus  is  available  on  a  regular  commer- 
ical  basis.  The  fixed  wing  equipment  operated  by  Renting 
Geophysics  in  Calgary  Canada  is  available  on  demand.  This 
system  is  described  by  Pitcher  et  al  (1930).  tfith  the 
exception  of  the  larger  coil  separation  and  of  necessity, 
larger  flight  altitude,  the  fixed  wing  system  operates  in  an 
identical  fashion  to  the  helicopter-towed  boom. 

il  Review  o f  Previous  Wor_k 

<*•2  E  lectr  omagnetlc  Measuremer.  ts 

There  are  no  known  attempts  at  measuring  the  thickness 
of  sea  ice  with  an  airborne  electromagnetic  induction 
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FREQUENCY  DOMAIN  SCHEMATIC 


COMMERICAL  RIGID  -  BOOM 

Airborne  EM  Systems 
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apparatus.  Keller  and  Frisehknecht  however  provide  an 
account  of  such  a  measurement  that  was  made  with  ground 
electr ©magnetic  equipment  over  thick  ice  in  the  Arctic.  The 
work  was  done  in  June  1959  (Keller  and  Frischknecht  1966, 
pp.  3U5-31*6)  in  an  area  where  the  ice  was  about  5  meters 
thick  and  the  water  had  a  conductivity  of  about  3.2  s/M. 
Eighty  soundings  were  made  using  a  coil  separation  of  the 
order  of  100m  and  a  primary  field  excitation  in  the  audio 
frequency  range.  Because  the  system  was  not  absolutely 
calibrated  readings  were  taken  at  a  number  of  frequencies 
and  a  curve  matching  procedure  was  used  to  interpret  the 
results.  On  the  average,  the  electromagnetic  determinations 
of  ice  thickness  appeared  to  be  accurate  to  about  10%. 

More  recently  Sinha  (1976)  conducted  a  limited  field 
study  to  evaluate  the  usefulness  of  a  light  weight,  small, 
high  frequency  portable  electromagnetic  prospecting  system 
for  the  determination  of  sea  ice  thickness.  Two  different 
pieces  of  commercial  equipment  were  used.  Both  had  a  coil 
separation  of  about  one  meter  but  differed  in  frequency  of 
operation  with  one  system  running  at  15  kHz  while  the  second 
one  operated  at  8  kHz.  Once  again  difficulties  were  encoun¬ 
tered  with  the  calibration  of  the  equipment  and  only  the 
vertical  gradient  of  the  system  response  could  be  properly 
measured.  The  ice  thickness  was  then  determined  from  this 
quantity.  Under  favourable  circumstances  accuracies  of  a 
few  centimeters  in  young,  thin,  (25-75m)  ice  could  be 
attained . 
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^.2  Sea  Ice  Properties 

Perhaps  the  most  reliable  set  of  data  on  the  electrical 
properties  of  sea  ice  was  obtained  by  McNeil  and  Hoekstra 
(1975)  when  they  attempted  to  measure  the  thickness  of  sea 
ice  by  measuring  the  wave-tilt  of  VLF  t r an sm i ss i on s . 
Although  this  work,  done  in  the  vicinity  of  the  Naval  Arctic 
Research  Laboratory  at  Point  Barrow,  Alaska  in  May  1971 
failed  to  prove  the  usefulness  of  the  wave-tilt  method  for 
the  desired  purpose,  it  did  result  in  a  fairly  complete  set 
of  data  on  the  electrical  properties  of  sea-ice.  Figure  6 
shows  a  summary  of  their  measurements.  The  data  clearly 
shows  that  young  first  year  ice  retains  much  of  its  salinity 
and  is  fairly  conductive  while  multi-year  ice  is  gradational 
in  resistivity.  As  expected,  multi-year  ice  shows  a  con¬ 
tinuous  progression  from  a  high  resistivity  of  about  10  K- 
Ohm-m  at  surface,  to  less  than  30  Ohm-m  at  the  sea  water 
surface. 

Based  on  the  Me Ne i 1 1-Ho e kstr a  data,  we  have  constructed 
some  sea  ice  models  which  were  then  used  to  evaluate  the 
feasibility  of  doing  a  test  of  airborne  electromagnetic  sea 
ice  thickness  determination.  Other  sea  ice  models  were  then 
used  to  determine  the  response  of  an  optimal  system  espe¬ 
cially  designed  for  sea  ice  thickness  mapping.  A  summary  of 
the  models  used  in  this  report  is  given  in  Table  2 . 
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TABLE  2 


Transitional  Models 
Multiyear  Ice  -  Thickness  T 


Model  1 

0 

0.2T 

p 

=  5009 n-M 

0.2T  - 

0.4T 

p 

=  1590o-m 

0.4T  - 

0.6T 

p 

=  200^-m 

0.6T  - 

0.8T 

p 

=  30«-m 

0.8T  - 

T 

p 

=  5«-m 

Model  2  0  - 

0 .251 

P 

=  100ft-M 

0.25T  - 

0.4QT 

P 

=  79^-m 

0.43T  - 

0.55T 

D 

=  45o-m 

0.55T  - 

0.70T 

P 

=  I5n-M 

0.70T  - 

0.85T 

P 

=  5ft-M 

0.85T  - 

T 

P 

=  ln-M 

Model  3 

0 

0.5T 

D  =  100^-M 

0.5T  - 

T 

o  =  10n-M 

Model  4 

0 

0.5T 

D 

=  100n-M 

0.5T  - 

0.75T 

P 

=  50^-m 

0.75T  - 

T 

P 

=  20^-m 

Model  5 

0 

0.5T 

P 

=  190n-M 

0.5T  - 

T 

P 

=  50^-m 

Model  6 

0 

T 

P 

=  33«-m 

r - — — 

i 

5 .  Determination  o f  Sea  Ice  Thickness  with  Airborne 

Electromagnetics 

5.  1  Survey  Principles  and  Choice  o f  System 

At  an  early  stage  of  this  research  project,  it  became 
evident  that  the  problem  of  mapping  ice  thickness  with  air¬ 
borne  electromagnetics  reduced  to  the  problem  of  determining 
precisely  the  vertical  distance  from  the  sensor  to  the  sur¬ 
face  of  the  conductive  sea  water.  Once  this  quantity  is 
known  then  one  simply  subtracts  from  it  the  vertical  dis¬ 
tance  from  the  sensor  to  the  ice-air  interface  as  determined 
by  a  precise  radar  or  laser  altimeter.  The  difference 
between  these  two  quantities,  i.e.  distance  between  sensor 
and  sea  surface  and  distance  between  sensor  and  ice  surface 
of  course  faithfully  represents  the  ice  thickness. 


It  is  unfortunately  impossible  to  accurately  establish 
the  sensor-conductor  distance  using  a  towed-bird  system 
without  knowing,  at  all  times,  the  position  of  the  bird  that 
contains  the  receiver.  Because  this  feature  is  not 
presently  available  on  any  towed-bird  system,  we  have 
automatically  rejected  the  use  of  a  time  domain  system  for 
the  purpose  of  determining  the  sea  ice  in  thickness  from  an 
aircraft.  Consequently,  all  subsequent  discussion  is  cen¬ 
tered  on  the  results  that  can  be  obtained  with  a  rigid-boom 
frequency  domain  system. 


5 . 2  F  easibility 

In  order  to  quickly  assess  the  feasibility  of  using  a 
helicopter  towed,  rigid-boom  electromagnetic  system  for  the 
determination  of  sea  ice  thickness,  we  computed  the  system 
response  for  a  device  with  a  co-axial  receiver  and  a 
transmitter  separated  by  a  distance  of  six  meters.  The 
results  are  presented  in  Figure  7  in  the  form  of  an  Argand 
diagram.  The  axes  are  conventionally  labeled  in  parts  per 
million  of  the  primary  field  at  the  receiver.  The  plot  is 
also  characterized  by  two  sets  of  parameters  namely  the 
vertical  distance  between  the  ir. strumentation  boom  and  the 
surface  of  the  conductor  and  "FC",  the  frequency- 
conductivity  product  for  the  conductor.  In  these  calcula¬ 
tions  it  was  assumed  that  the  sea  ice  if  any  was  infinitely 
resistive.  As  presented,  the  data  is  valid  for  any  6m  sys¬ 
tem.  For  example  if  we  use  a  hypothetical  6m  -  1000  hz  sys¬ 
tem  and  obtain  a  reading  of  -90 ppm  quadrature  and  -280ppm  in 
phase  (this  reading  is  located  at  point  "A"  on  the  chart 
shown  in  Figure  7)  then  we  can  deduce  that  the  sea  water 
conductivity  is  **  S/m  and  that  the  system  is  at  about  3'. 25m 
above  the  surface  of  the  water.  If  in  addition  to  this  we 
had  an  altimeter  reading  to  show  a  distance  of  30m  to  the 
ice  surface  then  an  ice  thickness  of  1.25m  can  be  inferred 
immediately . 

The  heavy  dashed  line  on  the  diagram  indicates  the  data 
locus  for  the  coaxial  Aerodat  system  operating  at  a  fre- 
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quenc.y  of  5U65  hz.  A  value  of  *»  S/m  is  assumed  for  the  sea 
water  conductivity.  At  this  point  it  is  worthwhile  to  note 
that  the  ratio  of  the  quadrature  to  the  in  phase  response  is 
fairly  constant.  Thus  the  phase  of  the  secondary  field 
varies  only  slightly  with  distance  to  the  sea  water  surface, 
the  amplitude  of  the  response  however,  is  very  strongly 
dependent  on  this  quantity. 

The  amplitude  of  the  secondary  field  for  the  Aerodat 
system  at  5^65  hz  above  •*  S/m  sea  water  is  plotted  as  a 
function  of  height  and/or  ice  thickness  in  Figure  8.  Xe 
note  a  gradient  of  about  30  ppm/m  of  ice  thickness.  Allow¬ 
ing  for  a  system  noise  of  1  ppm  one  can  immediately  conclude 
that  under  ideal  conditions,  this  system  could  detect  ice 
thickness  variations  of  as  little  as  3  cm! 

In  order  to  estimate  the  magnitude  of  error  that  might 
be  encountered  in  using  the  graph  shown  in  Figure  8,  which 
was  computed  under  the  assumption  of  infinitely  resistive 
ice,  to  interpret  data  that  is  related  to  ice  of  finite  con¬ 
ductivity  wa  have  computed  theoretical  data  for  ice  models  1 
and  6.  Model  6  represented  young  ice  with  an  average  resis¬ 
tivity  of  33  Oh.m-m,  Model  1  represented  multi-year  ice  with 
a  gradational  resistivity  profile.  The  theoretical  data  for 
these  models  were  then  "interpreted"  using  the  graph  of  Fig¬ 
ure  8.  The  results  of  this  operation  are  tabulated  ir.  Table 
3  and  clearly  indicate  that  finite  ice  resistivity  is  a  very 
minimal  source  of  error.  To  all  intents  and  purposes  sea 
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FIGURE  3  SYSTEM  RESPONSE 
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Errors  in  Ice  Thickness  Estimates 


ice  is  transparent  to  audio  frequency  electromagnetic 
energy . 

It  is  more  costly,  from  the  view  point  of  accuracy,  to 
ignore  the  effects  of  the  sea  water  conductivity  and  its 
depth.  The  simple  method  of  data  interpretation  that  relies 
on  the  ideal  data  shown  in  Figure  8  will  at  worst  result  in 
a  1 5 X  error  in  the  estimated  ice  thickness  for  a  50i  varia¬ 
tion  in  sea  water  conductivity.  Thus  if  the  sea  water  con¬ 
ductivity  is  3  S/m,  instead  of  the  supposed  **  S/m,  a  2  m 
layer  of  ice  will  appear  to  be  only  175cm  thick.  The 
effects  of  changes  in  sea  water  conductivity  however  can  be 
accounted  for  by  using  a  more  complicated  system  of  data 
interpretation  based  on  charts  such  as  that  shown  in  Figure 
7  which  allows  for  estimates  of  both  sea  ice  thickness  and 
water  conductivity. 

Neglecting  the  fact  that  the  sea  ice  can  be  floating  on 
an  inadequately  thick  layer  of  water  can  also  cause  very 
great  errors  in  the  estimates  of  ice  thickness.  In  fact, 
good  ice  thickness  estimates  can  only  be  obtained  where  the 
sea  water  is  at  least  one  skin  depth  (or  about  3.5  m  at  5^65 
Hz)  deep.  In  areas  where  the  water  is  shallower  the  effects 
of  bottom  conductivity  would  be  difficult  to  separate  out 
using  only  a  single  frequency  measurement. 

We  have  also  examined  the  possibility  of  using  the 
fixed  wing  Kenting/Seintrex  Tridem  system  that  operates  at 
8000  hz  and  normally  flies  at  about  50m  above  the  surface. 
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In  a  manner  analogous  to  Figure  8,  the  amplitude  response 


for  this  system  is  shown  in  Figure  9.  Here  we  note  a  gra¬ 
dient  of  about  250  ppm/m  of  sea  ice.  The  system  noise  how¬ 
ever  is  about  25  ppm  and  thus  its  resolution  is  only  about 
♦  10  cm  under  the  best  of  circumstances. 

5 . 3  A_£  Op  t  i  ma  1  AEM  System  for  Happing  the  Thickness  of 

Sea  Ice 

The  definition  of  an  optiml  airborne  electromagnetic 
system  for  mapping  the  thickness  of  sea  ice  can  best  be  done 
in  terms  of  a  number  of  specific  criteria.  We  must  have  a 
system  that  is; 

(  1  )  Simple  to  install 

(2)  Simple  to  operate 

(3)  Yields  data  that  can  be  interpreted  in  real  time. 

(M  Demonstrates  adequate  (£  10  cm)  sensitivity  to  ice 

thickness  variations. 

(5)  Is  relatively  insensitive  to  variations  in  other  param¬ 
eters  such  as  : 

a.  ice  resistivity 

b.  water  resistivity 

c.  water  depth 

(6)  Can  be  installed  on  a  fixed  wing  aircraft  to  allow  for 
missions  of  long  duration. 
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FIGURE  9  SYSTEM  RESPONSE 
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A  system  that  respons  to  all  these  criteria  could  be 
installed  on  a  slow  flying  modern  amphibian  aircraft  such  as 
the  Canadair  CL  I1***.  in  order  to  take  advantage  of  existing 
technology  we  suggest  a  horizontal  axis  transmitter.  Sim¬ 
plicity  of  operation  and  good  sensitivity  to  ice  thickness 
will  be  obtained  with  a  vertical  axis  receiver  mounted  in  a 
tail  stinger.  A  reasonable  transmitter-receiver  separation 
for  this  type  of  system  would  be  20  m.  It  could  be  flown  at 
50  m  above  the  ice  surface. 

The  frequency  of  operation  must  be  chosen  in  a  way  that 
renders  the  system  very  sensitive  to  the  position  of  the  sea 
water-ice  interface  but  not  to  the  electrical  properties  of 
either  the  sea  ice  or  the  sea  water.  In  particular  ore 
would  like  to  use  a  frequency  for  which  both  young  or 
multi-year  sea  ice  will  be  transparent  while  at  the  same 
frequency,  even  moderately  saline  sea  water  would  constitute 
a  perfect  reflector.  After  some  detailed  inspection  of 
these  criteria  30,000  Hz  appears  to  be  a  frequency  that 
satisfies  these  requirements. 

Under  these  assumptions,  we  can  readily  use  electromag¬ 
netic  image  theory  (Keller  and  Frischknecht  1966,  pp.  327- 
328)  to  predict  the  system  response.  This  theory  assigns  an 
"infinite"  conductivity  to  the  conducting  medium  so  that  the 
conductor  insulator  interface  (sea  water-sea  ice  in  our 
case)  forms  a  perfect  reflector  for  the  source.  In  spite  of 
this  approximation  this  theory  can  give  useful  results  in  a 


n umber  of  situations  (Duckworth,  197M 


In  our  case,  as  will  be  demonstrated  by  comparing  the 
exact  numerical  data  and  approximate  results  obtained  from 
image  theory,  this  theory  can  be  used  to  quite  accurately 
predict  the  optimal  system  response  if  the  computation  is 
done  as  fol lows  : 

Let  P  =  the  30  -<hz  in  phase  component  of  the  secondary 
field  . 

Let  Q  =  the  associated  quadrature  component 
Let  A  =  P  ♦  Q  (  1  ) 

denote  the  "image"  anomaly. 

Then,  as  suggested  by  the  relations  given  by  Keller  and 
Frischknecht  (1966,  p.  3U0),  "A"  as  defined  above  consti¬ 
tutes  an  excellent  approximation  to  the  ideal  "image  theory" 
system  response. 

For  the  optimal  system  in  question,  the  guiding  rela¬ 
tions  are: 

h  =  a  *  t  ( 2  ) 

h  =  system  distance  to  water  surface 
a  =  system  altitude  above  ice 
t  =  ice  thickness. 

If  "1"  denotes  the  system  coil  separation  then  the  system 
response  is  given  by: 


A  =  3tl/h)4  /[u+l^/h*  ]  ^  (3) 


The  geometrical  quantities  used  in  this  expression  are 
shown  in  Figure  10.  Figure  11  shows  the  system  response  as 
a  function  of  ice  thickness.  Allowing  for  the  fact  that  the 
ice  thickness  is  much  smaller  than  the  system  altitude  and 
that  the  coil  separation  is  also  small  in  comparison  with 
the  system's  distance  to  the  water  surface  it  is  possible 
to  immediately  derive  a  simple  approximate  relationship  for 
the  system  response.  This  turns  out  to  be: 

A  =  t  3/ 32  ]  1 1/ a  ]  ^"  [1-*»t/50]  (M 

For  the  suggested  optimal  system  (20m  coil  separation,  50m 
altitude)  equation  *»  simplifies  further  to; 

A  tr  2i*00(  1-i*t/50]  ppm  (5) 

This  expression  is  only  accurate  to  about  10t  when  com¬ 
pared  with  the  corresponding  values  obtained  from  equation 
3.  It  does  however  correctly  indicate  a  gradient  of  about 
200  ppm  per  meter  of  sea  ice.  Thus  if  we  are  to  resolve 
about  5  cm  of  ice  the  optimal  system  need  only  have  a  noise 
level  of  lOppm . 

We  next  examine  the  effects  of  sea  ice  conductivity  on 
the  correct  estimate  of  its  thickness  using  image  theory. 
This  is  done  in  Table  •*  which  illustrates  this  problem  for  a 
2  meter  thick  slab.  The  accuracy  of  the  image  method  is 
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OPTIMAL  SYSTEM  RESPONSE 


surprisingly  good  showing  the  worst  case  error  of  H»  ctn  for 
young  ice. 

The  effects  of  sea  water  conductivity  on  ice  thickness 
estimates  as  seen  for  a  2  meter  layer  of  33  ohm-ra  ice  are 
examined  in  Table  5.  Again,  the  proper  use  of  image  theory 
renders  this  effect  minimal  with  the  worst  error  amounting 
to  but  1**  cm  or  7%.  The  effects  of  sea  water  aepth  on  ice 
thickness  estimates  were  not  examined  for  the  optimal  sys¬ 
tem.  It  is  believed  however,  that  as  previously  mentioned 
the  water  must  be  at  least  one  skin  depth  deep  in  order  to 
correctly  estimate  ice  thickness.  At  the  high  operating 
frequency  of  the  optimal  system  this  requirement  is  trivial 
as  it  amounts  to  1.5m  of  sea  water. 

6 .  Conclusions  and  Recommendations 

On  the  basis  of  the  foregoing  discussion  one  can 
readily  conclude  that  it  is  quite  feasible  to  measure  the 
thickness  of  sea  ice  in  the  polar  regions  using  an  airborne 
e 1 ec tr ornagn e t i c  system.  In  order  to  avoid  the  errors  related 
to  receiver  motion  the  measurements  must  be  made  in  the  fre¬ 
quency  domain  usr.g  a  rigidly  coupled  apparatus.  In  princi¬ 
ple  either  a  helicopter  towed  boom  equipment  or  a  fixed  wing 
apparatus  can  be  used  for  the  purpose.  The  former  is  pre¬ 
ferred  for  small  scale  surveys  because  of  its  "button-on" 
adaptability  to  any  number  of  carrier  vehicles.  For  '  or. g 
range  large  scale  surveys  the  fixed  wing  version  is  prefer¬ 
able. 


31 


Influence  of  Seawater  Conductivity 


There  is  no  doubt  that  conventional  survey  equipment 
available  now  for  contract  work  would  give  valuable  results. 
It  car.  be  used  to  verify  the  theoretical  ideas  presented 
here  and  to  assure  us  that  no  possible  sources  of  error  have 
beer,  ig  cored.  There  is  little  to  choose  between  the  various 
commercial  systems  listed  in  Table  1. 

Finally,  we  conclude  that  the  simple  electr omagr.etic 
image  theory  can  be  used  to  predict  the  response  of  an 
optimal  system  for  mapping  the  thickness  of  sea  ice.  A  sys¬ 
tem  that  uses  this  technique  for  data  interpretation  car. 
readily  translate  the  observed  secondary  field  amplitudes  in 
terms  of  an  underlying  ice  thickness. 

Because  the  determination  of  sea  ice  thickness  with  an 
airborne  electromagnetic  device  is  entirely  feasible  we 
strongly  recommend  that  the  concepts  outlined  in  this  report 
be  subjected  to  a  comprehensive  field  test.  This  test  sould 
be  carried  under  operating  conditions  similar  to  those 
envisaged  for  future  practical  application  of  this  tech¬ 
nique.  We  suggest  the  use  of  a  helicopter  towed  system  that 
will  be  equipped  with  a  precise  altimeter  located  or.  the 
same  boom  as  the  electromagnetic  sensors.  Ground  verifica¬ 
tion  of  the  airborne  thickness  determinations  can  be  done  by 
a  number  of  drilled  holes.  Advantage  should  be  taken  of  the 
drilling  program  to  carry  out  a  limited  amount  of  measure¬ 
ments  on  the  electrical  properties  of  sea  ice. 


33 


RE  FERFNCES 


Barringer,  A.  R,  1962,  The  INPUT  electrical  pulse  prospecting 
system;  Min.  Cong.  J.,  vol.  **8,  pp.  **9-52. 

Becker,  A.  1969,  Simulation  of  time-domain,  airborne,  electro¬ 
magnetic  system  response.  Geophysics,  vol.  31*.  PP«  739-752. 

Becker,  A.  1979,  Airborne  electromagnetic  methods.  Geophysics 
and  Geochemistry  in  the  Search  for  Metallic  Ores;  Peter  J.  Hood, 
editor;  Geological  Survey  of  Canada,  Economic  Geology  Report 
31,  PP-  33-M. 

DeMoully,  G.T.  and  Becker,  A.,  1932,  Automated  interpretation  of 

airborne  electromagnetic  data.  Abstract.  Technical  Program, 
52nd  Annual  Mtg.  S.E.G.  pp.  382-381*. 

Duckworth,  K.  1970,  Electromagnetic  depth  sounding  applied  to 
mining  problems.  Geop.  vol.  35,  pp.  1086-1098. 

Fraser,  D.C.  1979,  The  multicoil  II  airborne  electromagnetic 
system.  Geophysics,  vol.  ****,  pp.  1367  —  1  39**. 

Keller,  G.V.  and  Fr i sc hkn ec ht ,  F.C.,  T  9  6  6 .  /Electrical 
methods  in  geophysical  prospecting,  Pergamon  Press,  N.Y. 

McNeill,  D.  and  Hoekstra,  P.  ,  1975,  In  situ  measurements  on 
the  conductivity  and  surface  impedance  of  sea  ice  at  VLF. 

Radio  Science,  vol.  8,  #1,  pp.  23-30. 

Morrison,  H.F.  and  Becker,  A.,  1982,  Analysis  of  airborne 

electromagnetic  systems  for  mapping  depth  of  sea  water. 

Final  report,  OMR  contract  ON  000 1 fc-82 -M -007 3 ,  57p. 

U . C . B . -Eng ineer i n g  Geoscience 

Pitcher,  D.H.,  Barlow,  R.B.,  and  Lewis,  M.,  1980,  Tridem 

airborne  conductivity  mapping  as  a  lignite  exploration 
method,  C1M  Bulletin,  May,  pp.  1-12. 

Sinha,  A.K.  (1976),  A  field  study  for  sea  ice  thickness 
determination  by  electromagnetic  means.  Geol.  Surv. 

Can.  Paper  76-1C,  pp.  225-228. 


34 


Distribution  List 


Department  of  the  Navy 
Asst  Deputy  Chief  of  Navy  Materials 
for  Laboratory  Management 
Rm  1062  Crystal  Plaza  Bldg  5 
Washington  DC  20360 

Department  of  the  Navy 
Asst  Secretary  of  the  Navy 
(Research  Engineering  &  System) 
Washington  DC  20350 

Project  Manager 
ASW  Systems  Project  (PM-4) 
Department  of  the  Navy 
Washington  DC  20360 

Department  of  the  Navy 
Chief  of  Naval  Material 
Washington  DC  20360 

Department  of  the  Navy 
Chief  of  Naval  Operations 
ATTN:  OP  951 
Washington  DC  20350 

Department  of  the  Navy 
Chief  of  Naval  Operations 
ATTN:  OP  952 
Washington  DC  20350 

Department  of  the  Navy 
Chief  of  Naval  Operations 
ATTN:  OP  987 
Washington  OC  20350 

Director 

Chief  of  Naval  Research 

ONR  Code  420 

ATTN:  Dr.  G.  B.  Morris 

Ocean  Science  &  Technology  Det 

NSTL,  MS  39529 

Director 

Defense  Technical  Info  Cen 
Cameron  Station 
Alexandria  V A  22314 


Commander 

DW  Taylor  Naval  Ship  RSD  Cen 
Bethesda  MD  20084 

Commanding  Officer 
Fleet  Numerical  Ocean  Cen 
Monterey  CA  93940 

Director 

Korean  Ocean  R&D  Inst 
ATTN:  K.  S.  Song,  Librarian 
P.  0.  Box  17  Yang  Jae 
Seoul  South  Korea 

Commander 

Naval  Air  Development  Center 
ATTN:  Dr.  A.  R.  Ochadlick,  Jr. 
Warminster  PA  18974 

Commander 

Naval  Air  Systems  Command 
Headquarters 

ATTN:  Code  340J  (Mr.  Barry  Dillon) 
Code  330G  (Dr.  Paul  F. 
Twitchell ) 

Washington  DC  20361 

Commanding  Officer 

Naval  Coastal  Systems  Center 

Panama  City  FL  32407 

ATTN:  Code  4130  (Dr.  G.  Kekelis) 

Commander 

Naval  Electronic  Sys  Com 
Headquarters 
Washington  DC  20360 

Commanding  Officer 
Naval  Environmental  Prediction 
Research  Facility 
Monterey  CA  93940 

Commander 

Naval  Facilities  Eng  Command 

Headquarters 

200  Stovall  St. 

Alexandria  VA  22332 


k 


Commanding  Officer 
Naval  Ocean  RAO  Activity 
ATTN:  Codes  UO/1L1 
Code  125 
Code  200 
Code  300 
Code  115 
Code  500 
NSTL  MS  39529 

Director 
Liaison  Office 
Naval  Ocean  RAD  Activity 
800  N.  Quincy  Street 
502  Ballston  Tower  #1 
Arlington  V A  22217 

Commander 

Naval  Ocean  Systems  Center 
San  Diego  CA  92152 

Commanding  Officer 
Naval  Oceanographic  Office 
ATTN:  Code  8000 
Code  8200 
Code  8400 
Code  8402 
NSTL  MS  39522 

Commander 

Naval  Oceanography  Command 
ATTN:  Larry  Reshew 
NSTL  MS  39522 

Superintendent 
Naval  Postgraduate  School 
ATTN:  Dr.  0.  Heinz,  Code  61  Hz 
Monterey  CA  93940 

Commanding  Officer 
Naval  Research  Laboratory 
Washington  DC  20375 

Commander 

Naval  Sea  System  Command 
Headquarters 

ATTN:  Daniel  Porter,  Code  63  R1 
Washington  DC  20362 

Commander 

Naval  Surface  Weapons  Center 
Dahlgren  VA  22448 


Commanding  Officer 
Naval  Underwater  Systems  Center 
ATTN:  New  London  Lab 
Newport  RI  02840 

Director 

New  Zealand  Oceano  Inst 
ATTN:  Library 
P.  0.  Box  12-346 
WELLINGTON  N.,  NEW  ZEALAND 

Director 

Office  of  Naval  Research 
Ocean  Science  A  Technology  Div 
NSTL  MS  39529 

Department  of  the  Navy 
Office  of  Naval  Research 
ATTN:  Code  102 
800  N.  Quincy  St. 

Arlington  VA  22217 

Commanding  Officer 
ONR  Branch  Office 
536  S  Clark  Street 
Chicago  IL  60605 

Commanding  Officer 
ONR  Branch  Office  LONDON 
Box  39 

FPO  New  York  09510 

Commanding  Officer 
ONR  Western  Regional  Ofs 
1030  E.  Green  Street 
Pasadena  CA  91106 

President 
Texas  AAM 

ATTN:  Dept  of  Ocean  Working  Collection 
College  Station  TX  77843 

Director 

University  of  California 
Scripps  Institute  of  Oceanography 
P.  0.  Box  6049 
San  Diego  C a  92106 

Director 

Woods  Hole  Oceanographic  Inst 
Woods  Hole  MA  02543 


36 


Chief  of  Naval  Material 
(MAT-0724) 

ATTN:  Capt.  John  Harlett 
800  N.  Quincy  St. 

Arlington,  VA  22217 

EUREKA  Resources  Associates,  Inc. 
ATTN:  R.  0.  Prindle 
2161  Shattuck  Avenue,  Suite  317 
Berkely,  CA  94704 

SRI  International 

Electromagnetics  Science  Laboratory 
ATTN:  A.  M.  Buberie 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

Rockwell  International 
Missile  System  Division 
ATTN:  S.  Stasenko 
3370  Miraloma  Avenue 
Anaheim,  CA  92803 

Geological  Survey  of  Canada 
Terrain  Geophysics  Subdivision 
ATTN:  L.  S.  Collett 
601  Booth  St. 

Ottawa,  Ontario 
Canada  K1A0E8 

Dobrocky  Seatech  Limited 
ATTN:  Dr.  J.  R.  Harper 
9865  West  Suanich  Rd. 

P.  0.  Box  6500 
Sidney,  British  Colombia 

Geophex,  Ltd. 

ATTN:  Dr.  I.  J.  Won 
319  Morrison  Avenue 
Raleigh,  NC  27608 

International  Exploration,  Inc. 
ATTN:  T.  L.  Jacobsen 
577  Sackettsford  Road 
Warminster,  PA  18974-1398 


Carson  Geoscience  Company 
ATTN:  A.  J.  Navazio 
32-H  Blooming  Glen  Road 
Perkasic,  PA  18944 

Airborne  Systems,  Inc. 

ATTN:  R.  L.  Summer felt 
10700  Walker  St. 

Cypress,  CA  90630 

Westinghouse  Electric  Corporation 

Defense  Group 

ATTN:  H.  L.  Hamrick 

98  Miracle  Strip  Parkway 

Ft.  Walton  Beach,  FL  32548 

EG&G  Washington  Analytical  Service 
Center,  Inc. 

Attn:  G.  L.  Bill 
2150  Fields  Rd. 

Rockville,  MD  20850 

Geotech,  Ltd. 

ATTN:  E.  B.  Morrison 
21-101  Amber  Street 
Markham,  Ontario 
Canada  L3R  3B2 

Barringer  Research  Ltd. 

ATTN:  C.  Cumming 
304  Carlingview  Drive 
Rexdale,  Ontario 
Canada  M9W  5G2 

Dynamics  Technology,  Inc. 

ATTN:  S.  R.  Borchandt 
22939  Hawthorne  Blvd. 

Torrance ,  CA  90505 

Chief  of  Naval  Research 
ONR 

ATTN:  Dr.  C.  A.  Luther,  Code  425AR 
Dr.  T.  Warfield,  Code  220 
Washington,  DC  20350 


37 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whtn  Palm  gnl.f.dj 

I  REPORT  DOCUMENTATION  PAGE 


^REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

NOR  DA  Technical  Note  261  yy  (? 

4.  TITLE  (mnd  Subtitle) 

ANALYSIS  OF  AIRBORNE  ELECTROMAGNETIC 
SYSTEMS  FOR  MAPPING  THICKNESS  OF  SEA  ICE 


NORDA  Technical  Note  261 


READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 


5  TYPE  OF  REPORT  6  PERIOD  COVERED 


r.  author^ j; 

A.  Becker 
H.  Morrison 
K.  Smits 

*•  PERFORMING  ORGANIZATION  NAME  AMD  ADDRESS  " 

Naval  Ocean  Research  and  Development  Activity 
Ocean  Science  Directorate 

NSTL,  Mississippi  39S29  _ _ 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


6.  PERFORMING  ORG.  REPORT  NUMBER 


ft.  CONTRACT  OR  GRANT  NUMBERS) 


10-  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


62759N _ 

12.  REPORT  DATE 

November  1983 

II  NUMBER  OF  PAGES 


Same _ 2fi _ 

14  MONITORING  AGENCY  NAME  4  AOORESSf//  dliiofnt  from  Comrol/lnfl  Oftlcm)  15  SECURITY  CLASS,  (of  lltlt  r.port; 

UNCLASSIFIED 

15.  declassification  downgrading 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Rmport) 

Approved  for  Public  Release 
Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (ol  (h«  mbotrmct  motored  In  Block  20.  II  dlllmrmnt  Iro m  Rmport) 


*».  SUPPLEMENTARY  NOTES 


1ft-  KEY  WOROS  (Continum  on  rmvmrmm  mldo  II  nmcmmmmry  an d  Idmntlfy  by  block  numbor) 

ice  thickness 

airborne  electromagnetic  (AEM) 
airborne  mapping 
geophysical  prospecting 

20.  ABSTRACT  fContfnu*  on  rmvmrmm  mldm  II  nmcmmmmry  end  Idmntlly  by  block  micnborj 

NORDA  Code  372  is  at  present  engaged  in  research  projects  for  extending  Air¬ 
borne  Electromagnetic  (AEM)  technology  from  geophysical  prospecting  applica¬ 
tions  to  Navy  requirements  such  as  airborne  shallow  water  bathymetry,  airborne 
ice  thickness  determination,  and  conductivity  profiling.  This  report  is  part  of 
our  research  program  addressing  AEM  ice  thickness  determination  and  will  demon¬ 
strate  that  the  Navy  requirements  for  rapid  airborne  methods  for  mapping  ice 
thickness,  in  real  time,  can  be  accomplished  with  Aem  techniques.  At  present 


00  I  JAM *71  1473  COITION  OF  I  NOV  ••  11  OMOLCTC 

S  'N  0102-  LF-  014-  6401 


UNCLASSIFIED 

ICCUNITY  CLASSIFICATION  OF  TMIC  ^A«l 


UNCLASSIFIED _ 

HtCUWTV  CLASSIFICATION  OF  THI»  NAO«  D**»  lnMn4  _ _ 

I  there  are  no  reliable  airborne  techniques  for  measuring  sea  ice,  although  impulse 
|  radar  has  yielded  promising  results  over  fresh-water  ice  but  not  over  sea-water 
ice. 


S/N  0102-  LF.  014-6601 

UNCLASSIFIED 


StCUNlTY  CLASSIFICATION  OF  THIS  FAOtfWMa  Oat*  AhmO 


